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I. INTRODUCTION 



The CPT theorem requires that the proton and antiproton Cp) lifetimes are equal. Searches for proton decay have 
yielded lower limits on the proton lifetime r p > O(10 32 ) yr Q. Hence, a search for p decay with a short lifetime 
(Tp << t p ) tests the CPT theorem. In this paper we describe a search for p decay at the Fermilab Antiproton 
Accumulator. 

CPT invariance is one of the most fundamental symmetries of modern physics. Over the years there have been 
a variety of searches Jl]-[| for CPT violation based on comparing particle and antiparticle masses, lifetimes, and 
magnetic moments. For example, the proton and antiproton masses have been shown to be equal with a precision 
of a few parts in 10 s ||, whilst the particle and antiparticle masses in the neutral kaon system have been shown to 
be equal to about one part in 10 18 A search for p decay complements these CPT tests by providing a search 

for CPT violation accompanied by a violation of baryon number. Indeed, since the antiproton is the only long-lived 
antiparticle that could in principle decay into other known particles without violating charge conservation, a search 
for p decay provides a unique test of the CPT theorem, and a unique test of the intrinsic stability of antimatter. 

The sensitivity of a p decay search to the presence of a CPT violating interaction has been characterized || by 
considering a dimension-6 CPT- violating operator with characteristic mass scale mj. Dimensional analysis then 
provides the estimate Tp ~ m\/m^ yielding: 



mx 



(4.3 x 10 9 GeV) x (rp/lO Myr) 1/4 . (1) 



Hence, a search for p decay with a lifetime of order 1 yr, for example, provides a test for CPT violation at a mass 
scale O(10 7 ) GeV/c 2 , which is well beyond the scale accessible at high energy colliders. 

The most stringent lower limit on Tp has been obtained || from a comparison of recent measurements of the cosmic 
ray p flux Q with predictions based on expectations for secondary production of antiprotons in the interstellar medium. 
The agreement between the observed and predicted rates implies that Tp is not small compared to T/7, where T is 
the p confinement time within the galaxy (~ 10 7 yr) and 7 is the Lorentz factor for the observed antiprotons. After 
taking into account the uncertainties on the relationship between the interstellar p flux and the flux observed at the 
Earth, at 90% C.L. the limit Tp > 8 x 10 5 yr has been reported S. This indirect limit is not valid if current models 
of p production, propagation, and interaction in the interstellar medium are seriously flawed. Indeed, it has been 
claimed || that, within the framework of the minimal supersymmetric extension to the standard model, if neutralino 
annihilations contribute significantly to the flux of cosmic ray antiprotons in the galactic halo, p lifetimes as low as 
10 5 years can accommodate the current cosmic ray data. 

Laboratory searches for p decay have, to date, provided less stringent limits on Tp. However, these limits do not 
suffer from large model dependent uncertainties. The most stringent published laboratory limit on inclusive p decay 
has been obtained from a measurement of the containment lifetime of ~ 1000 antiprotons stored in an ion trap, 
yielding T-p > 3.4 months The sensitivity of laboratory p decay searches can be dramatically improved by looking 
for explicit p decay modes at an antiproton storage ring. Angular momentum conservation requires that a decaying p 
would produce a fermion (electron, muon, or neutrino) in the final state. Hence, there are only a limited number of 
potential decay modes. A search for explicit p decay modes with an electron in the final state was made by the APEX 
test experiment (T861) at the Fermilab Antiproton Accumulator. The T861 search yielded the 95% C.L. limits [||: 
Tp/B(p -> e"7) > 1848 yr, Tp/B{p -> e-n°) > 554 yr, Tp/B{p -> e'-q) > 171 yr, T ¥ /B(p -> e~ K%) > 29 yr, and 
T ¥ /B(p^e-K° L )> 9 yr. 

Following the T861 test experiment, the APEX detector |||l0| was designed to enable a more sensitive search for 
p decay at the Fermilab Antiproton Accumulator. Based on APEX data, we have recently reported limits on Tp/B 
from the first search for antiproton decay modes with a muon in the final state [ pi||l3j ]. In this paper we give a more 
detailed description of our limits on p decay modes with a muon in the final state, and describe new results from a 
search for p decay modes with an electron in the final state. 

The paper is organized as follows. Sec. II gives a brief description of the APEX detector. The APEX triggers and 
data samples are described in Sec. III. The search for decay modes with a muon in the final state is described in 
Sec. IV., and the electron decay mode search is presented in Sec. V. Conclusions are summarized in Sec. VI. 



II. THE APEX DETECTOR 



The APEX detector was located in a 15.9 m straight section of the Fermilab Antiproton Accumulator ring, which 
has a circumference of 474 m. The detector was designed to identify p decays within a short 3.7 m long section of 
the ring, and was optimized to detect a single energetic charged track (electron or muon), originating from the beam, 
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and accompanied by one or more neutral pions or photons. A schematic of the detector is shown in Fig. In the 
following we give a brief description of the main detector components relevant to the analysis described in this paper. 
A more detailed description can be found in Ref . (To) . The APEX detector consisted of the following: 

(i) A 3.7 m long evacuated decay tank operated at 10" 11 Torr. The vacuum was maintained by a roughing pump, an 

ion diffusion pump, and twelve titanium sublimation filaments. The tank was designed to minimize the residual 
gas within the fiducial volume of the experiment, and hence minimize backgrounds from beam-gas interactions. 
The downstream section of the tank consists of a 96 cm diameter cylindrical shell supporting a 1.2 mm thick 
(0.073 radiation lengths) stainless steel window. 

(ii) A movable tungsten wire target at the upstream end of the tank (z = where, in the APEX coordinate system, 

z is measured in the direction of the antiproton beam). The target could be inserted into the beam halo to 
produce a localized source of particles for aligning and calibrating the detector. 

(iii) An upstream veto system consisting of horizontal and vertical scintillation counters arranged around the 10 cm 
diameter beam pipe immediately upstream of the tank and target. The counters covered a 1 x 1 m 2 area normal 
to the beam direction, and were used to reject tracks from upstream interactions. 

(iv) Three planes of horizontal and three planes of vertical scintillation counters downstream of the tank. Each 
plane consisted of two 50 x 100 x 1.27 cm 3 counters. The last planes of horizontal and vertical counters were 
downstream of a 2.3 radiation length lead wall, providing a preradiator to aid in identifying electrons and 
photons. The remaining counter planes were upstream of the lead, and provided pulse height information used 
to count the number of charged particles in an event (dE/dx counters). The measured pulse height distribution 
arising from the passage of minimum ionizing particles (MIPs) through the dE/dx counters is well described by 
a GEANT Monte Carlo simulation in which the yield is 37 photoelectrons per MIP. The predicted efficiency for 
a minimum ionizing particle to produce a counter signal in excess of a 0.5 MIP threshold is 97%. 

(v) A lead-scintillator sampling electromagnetic calorimeter [|ll] constructed from 144 rectangular 10 x 10 cm 2 

modules that are 14.7 radiation lengths deep. The modules were arranged in a 13 x 13 array with 6 modules 
removed from each of the four corners, and the central module removed to allow passage of the beam pipe 
(Fig. ||). Each module was read out with one 5.1 cm diameter photomultiplier. An LED pulser system enabled 
the stability of the photomultiplier gains for each module to be monitored. The calorimeter was calibrated by 
measuring the response to minimum ionizing tracks and reconstructing 7r° — > 77 and r\ — > 77 mass peaks using 
data samples recorded with the calibration target inserted into the beam halo. The measured mass peaks have 
fractional rms widths given by cr m /m ~ 0.25 for cluster pairs in the energy range of interest for the experiment. 
This mass resolution is dominated by the energy resolution of the calorimeter. 

(vi) A tail catcher (TC) downstream of the calorimeter consisting of a 20 cm deep lead wall followed by planes of 
horizontal and vertical scintillation counters. 

(vii) A limited-acceptance muon telescope (MT), 10 nuclear interaction lengths deep, located downstream of the TC, 
and aligned to point towards the center of the decay tank. The MT consists of a sandwich of five iron plates and 
five 30 x 30 cm 2 scintillation counters. The MT was used to identify charged particles that penetrate through 
several interaction lengths of material (muon candidates). 

(viii) A tracking system consisting of three planes of horizontal and three planes of vertical 2 mm diameter scintillating 
fibers downstream of the tank and upstream of the preradiator lead. The fibers were measured to yield about 
9 photoelectrons for a traversing minimum ionizing particle, and were read out using 12-stage Hamamatsu 
multianode photomultipliers. These detectors provided three space points along the track trajectory with 
typical residuals of 620 /im in the directions transverse to the beam direction. This resolution enabled the origin 
of tracks emerging from the decay tank to be reconstructed with an rms precision given by a z = 12 cm. The 
average single-hit efficiency for a minimum ionizing particle passing through a tracking plane was 85%. The 
measured track reconstruction efficiency, for a sample of events that are consistent with having one minimum 
ionizing particle passing through the dE/dx planes and the MT, is (90 ± 7)%. 

III. DATA TAKING 

The APEX experiment took data during the April-July 1995 Fermilab Tevatron collider running period at times 
when there were typically 10 12 antiprotons circulating in the Accumulator ring with a central p momentum of 8.90 ± 
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0.01 GeV/c (7 = 9.54±0.01). APEX ran parasitically to the Tevatron collider program. Due to increased background 
interactions when antiprotons were being stacked in the Accumulator, data were only taken during periods when 
stacking was not taking place. Such times occurred during the short period (typically 90 minutes) before daily 
injection of antiprotons into the Tevatron, or at times when the collider complex was not running due to maintenance 
or repairs. A measure of the sensitivity of the APEX data sample is given by: 



where Np-(t) is the number of circulating antiprotons at time t, the integral is over the live-time of the experiment, 
and the uncertainty arises from the precision with which the time dependent beam current was recorded. Hence, if 
Tp = 3.31 x 10 9 years then on average 0.63 antiproton decays would have occurred somewhere within the Accumulator 
during the live-time of the experiment. 

Energetic particles passing through the detector during Accumulator operation predominantly arise from interac- 
tions of the p beam with the residual gas in the decay tank or with material surrounding the beam. To suppress these 
backgrounds, and select candidate p — > [i~X and p — > erX decays, signals from the calorimeter and the scintillation 
counters were used to form muon and electron triggers. 



To search for muonic p decays, data were recorded with a muon trigger that required a coincidence between at 
least two of the five MT scintillation counters. These triggers were eliminated if they were in coincidence with a 
signal in one or more of the upstream veto counters indicating the presence of an interaction upstream of the decay 
tank. This loose trigger resulted in 1.2 xlO 6 events being recorded with a typical trigger rate of ~ 1 Hz. These 
events predominantly arise from background interactions in which the coincident MT counter signals are produced by 
traversing muons coming from charged pion decays, or by hadronic showers not contained in the calorimeter and TC. 



To search for electronic p decays, data were recorded with three complementary triggers: (i) An Et, a i ance trigger 
designed to select events consistent with a two-body p decay that results in the deposition of energy in two diagonally 
opposite quadrants of the calorimeter. The calorimeter cell groupings into four trigger quadrants are shown in Fig. |2j. 
The trigger required the summed signals in each of two diagonally opposite quadrants to exceed a threshold that was 
set equal to about 2 GeV; (ii) An E to tai trigger designed to select p decays in which all of the decay products deposit 
approximately all of their energy within the calorimeter. This would result in a total calorimeter energy deposition 
Etotai close to the beam energy. The trigger therefore required the summed signals from all of the calorimeter cells 
to exceed a threshold that was set equal to about 7 GeV; and (iii) An electron trigger which required signals in 
the scintillators of the first two dE/dx planes consistent with the passage of a MIP, a signal in the corresponding 
preradiator scintillators exceeding about 1.5 MIPs, and a signal in the corresponding calorimeter quadrant exceeding 
a summed energy threshold set equal to about 5 GeV. 

All three triggers were combined with a signal from the upstream veto counters, used in anti-coincidence to suppress 
events from upstream interactions. The typical data taking rates from background events satisfying the E oa i ance , 
Etotai, and electron trigger requirements were respectively 3 Hz, 9 Hz, and 8 Hz when there were 10 12 antiprotons 
circulating in the accumulator. The E oa i ance , Etotai, and electron triggers resulted in a data sample of 37.8 million 
events recorded to tape for further analysis. 



A number of additional triggers were used, either during normal data taking or during special running periods, to 
collect data for monitoring and calibrating the detector. Three of these triggers provided calibration data samples that 
were specially important for the analysis: (1) A very low threshold rate limited calorimeter trigger [Etotai > 250 MeV) 
that recorded events with minimal bias continuously during normal running. Dedicated minimal bias data samples 
were also taken during extended periods when there was a low intensity p current (corresponding to ~ 10 11 stored 
antiprotons) in the Accumulator; (2) A traversing charged track trigger that required a coincidence between the 
upstream veto counters and the TC counters. This trigger provided a sample of MIPs passing through the calorimeter; 




(2) 



A. Muon trigger and data sample 



B. Electron trigger and data samples 



C. Calibration triggers 
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and (3) Low threshold (E to tai > 5 GeV) and high threshold (E to tai > 7 GeV) calorimeter triggers used during special 
periods during which the calibration target was inserted into the beam halo. This provided tracks coming from a 
known origin which facilitated the relative alignment of the tracking planes. These data also provided a sample 
of localized interactions which enabled the reconstruction of the ir° — > 77 and 77 — > 77 peaks used to define the 
calorimeter calibration. 



IV. SEARCH FOR ANTIPROTON DECAY WITH A FINAL STATE MUON 



We begin with the sample of 1.2 xlO 6 events recorded with the muon trigger. Offline, after final calibration of 
the scintillation counters, the upstream veto counter requirement was re-applied using a more stringent threshold. 
This reduced the data sample to 1.1 xlO 6 events. The remaining backgrounds can be further suppressed by requiring 
that (a) the scintillation counter and scintillating fiber tracker signals are consistent with the presence of a minimum 
ionizing charged track that points back to the beam trajectory within the decay tank and that points forward to a 
pattern of minimum ionizing energy depositions within the MT counters, and (b) an event topology and kinematics 
that is consistent with the decay p — » \i~ + X. In the following, we consider first the simplest decay channel p — > fi~j, 
and then describe the search for more complicated decay modes. 



A. p — > /i 7 event selection 



To select events consistent with the decay p — > /i _ 7 we begin by requiring that: (i) at least four of the five 
MT counters be above threshold (4.2 x 10 4 events), and (ii) there be one and only one reconstructed track that 
extrapolates to the MT counters within the expected uncertainty due to multiple scattering, and also extrapolates 
back to the beam orbit with a point of closest approach within the fiducial volume of the decay tank (0 < z < 375 cm) 
and with an impact parameter |l(| less than 1 cm (416 events). Since a p — ► fi~j decay would result in events 
in which an energetic photon is produced coplanar with the muon (i.e. traveling within the plane defined by the 
muon and the incoming beam direction), we next require that (iii) there be at least one neutral cluster |l0| ] in 
the calorimeter (209 events) that is coplanar with the muon candidate (±5°). The data sample is reduced to 14 
events. The observed neutral cluster energy distribution for these events is compared in Fig. |3| with the predicted 
distribution obtained using the GEANT |l4j simulation described below, and corresponding to p — * /x~7 decay with 
a lifetime Tp/B(p — > (J,~j) = 5000 yr. The observed distribution peaks at low cluster energies, with a tail extending 
to approximately 2 GeV. In contrast, the predicted distribution for p — » /j~7 decays peaks at about 3.5 GeV, with 
only 4.4% of the simulated events having cluster energies less than 2 GeV. We conclude that there is no evidence for 
a signal. 



B. Limit on p — > p, 7 decays 



To set a limit on Tp/B(p — > /Lt _ 7) the remaining background events are eliminated by requiring that the neutral 
cluster energy exceed a minimum value E m i n , where we choose E m i n = 3 GeV. We note that although the choice of 
E m in is somewhat arbitrary, the final limit that we obtain on the decay p — + /_i~7 is insensitive to small changes in 
E m i n . After setting E m i n = 3 GeV, the resulting limit on Tp/B(p — ► /U~7) is given in years by: 

T p/B(P M _ 7) > - ln(1 ^ max /eSy (3) 

where N max = 2.3 is the 90% C.L. upper limit on the observation of N — events, and e is the fraction of decays 
taking place uniformly around the accumulator ring that would pass the trigger and event selection requirements. If 
{N max /eS) << 1, to a good approximation : 

T- p /B{p^ iri) > - f N ¥ (t)dt = (3.31 ±0.03) x 10 9 . (4) 

7 1* max J ^ max 

To take account of a e , the systematic uncertainty on e, we use the prescription given in Ref. p5| , giving at 90% C.L. 

N ma x = l^max X (1 + fJ'max "r/^)' ("-*) 

where a r = c e /e, and fj, m ax is the upper limit corresponding to the observation of N events (fJ, max = 2.3 for N = 0). 
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The GEANT Monte Carlo program has been used to simulate the detector response and calculate e. The detector 
simulation includes a full description of the detector geometry, and correctly describes the calorimeter, tracker, and 
scintillation counter responses (dE/dx, MT) measured using calibration data, together with the measured performance 
of the muon trigger. Further details can be found in Rcf. |lC(]. The efficiency e was obtained by generating 10 5 p — > /i~7 
decays uniformly along the beam orbit within the decay tank, and applying the trigger and data selection requirements 
to the simulated events. The geometrical efficiency of the detector and trigger is negligible for decays occurring outside 
of the tank. We obtain e = (3.7 ± 0.9) x 10 -5 , where the uncertainty on e arises from the systematic uncertainties on 
the muon trigger counter threshold calibration, the calorimeter energy scale, and the track reconstruction efficiency. 
The contributions to ay arising from the uncertainties on the muon trigger counter threshold calibration and the 
calorimeter energy scale are, respectively, ±18% and ±15%. These uncertainties have been evaluated by analyzing 
GEANT Monte Carlo p — > pT^j samples in which the simulated trigger calibration and the calorimeter energy scale 
have been changed by ±1 standard deviation. The overall systematic uncertainty on e was calculated by adding 
these contributions in quadrature with the uncertainty on the track reconstruction efficiency (±7%). Inserting the 
calculated e and a e into Eqs. 3 and 5 we obtain the result @ Tp/B(p -> ^-7) > 5 x 10 4 yr (90% C.L.). 

Finally, we note that although our limit is insensitive to the exact value of E m i n , the dependence of the final result 
on the position of this last cut is somewhat unsatisfying. Hence, to verify our result, we have also used the more 
complicated statistical method described in Ref. (r|. In this alternative statistical analysis all of the events in the 
final (E m in) distribution are retained, and use is made of how "signal-like" the position of each event is in the final 
distribution. Using this method we obtain a limit on Tp/B(p — > /i~7) that is consistent with our quoted result. 

C. Search for p — » fi~ir°, p — » pL~ri, and p — > /X - 77 decays 

Now consider the two-body decays p — > /i~7r° and p — > fjTT] (it, rj — ► 77), and the three-body decay p — > /i~77- 
These decays would result in events with one or two neutral clusters observed in the calorimeter, where the one-cluster 
events occur when the two photon-showers are not spatially resolved in the calorimeter or when one of the photons is 
outside of the calorimeter acceptance. In addition, the 77 may also decay into more complicated final states producing 
further clusters in the calorimeter. To optimize the search for /i~7r° and pTr\ final states we divide the 209 events 
described previously that have a muon candidate plus one or more neutral clusters into two subsamples, namely a 
one-cluster sample containing 104 events, and a multi-cluster sample containing 105 events. The position of the 
cluster in the one-cluster sample is required to be coplanar with the muon (±10°). The position of the cluster-pair 
formed by the energy-weighted centroid of the two highest energy clusters in the multi-cluster events is also required 
to be coplanar with the muon (±10°). These requirements reduce the samples to 13 one-cluster and 10 multi-cluster 
events. Note that we expect the multi-cluster events to contain 2 photons if they are genuine p — > /U _ 7T , fx" rj, or 
/i~77 decays. In most cases at least one of the two photons would be expected to convert in the preradiator lead. We 
therefore require that the preradiator counters in the quadrant opposite the MT have a signal exceeding a threshold 
of 0.5 MIPs. Only 3 multi-cluster events satisfy this requirement. Hence we are left with 16 events with > 1 cluster 
(fi + > 1 cluster events) for further analysis. 

Under the hypothesis that the observed (n'y) and (^77) systems arise from the decay of a beam particle, the mass 
of the parent particle can be computed from the measured muon direction and the directions and energies of the 
neutral clusters, using the constraint that the vector sum of the momentum components of the daughter particles 
transverse to the beam direction is zero. The resulting mass distribution for the remaining 16 (p + > 1 cluster) 
events is compared in Fig. ^ with predictions from the GEANT simulation for the decays (a) p — » H~tt°, (b) p — » fJ,~r), 
and (c) p — > /i~77. The observed mass distribution peaks at low masses with a tail extending to approximately 
0.7 GeV/c 2 . In contrast to this, the simulated signal distributions peak at the p mass, with only 17% (23%) [13%] of 
the simulated fi~ir° (pTr)) [/i _ 77] decays resulting in reconstructed masses less than 0.7 GeV/c 2 . We conclude that 
there is no evidence for a signal. 

To set limits on Tp/B we eliminate the remaining background by requiring that the reconstructed mass exceed a 
minimum value m m i n , and choose m m i n = 0.75 GeV/c 2 . We note that although the choice of m m i„ is somewhat 
arbitrary, the final limits that we obtain on the decays p — > p~Tr°, ^i~r], and p~^fj are insensitive to small changes 
in m m i n . The calculated overall efficiencies e for these decays to pass our trigger and analysis requirements are 
(3.6 ±0.9) x 10~ 5 for the p~tt° mode, (6.1 ±1.9) x 10~ 6 for the pTr\ mode, and (1.8 ±0.5) x 10~ 5 for the ^"77 mode. 
Substituting the values into Eqs. 3 and 5 yields the limits Tp/B(p — > h~tt°) > 5 x 10 4 yr, Tp/B(p — > (J,~r]) > 
8 x 10 3 yr, and T ¥ /B(p -> ^"77) > 2 x 10 4 yr at 90% C.L. 

Finally, to verify our result, we have removed the last cut on our final (mass) distribution, and used the more 
complicated statistical method described in Ref. p6|| . The resulting limits on the /x _ 7r°, /i - ?7, and /^ _ 77 decay modes 
are consistent with our quoted results. 
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D. Other p — ► fi + X decays 



Now consider the other possible two-body muonic decay modes of the antiproton, namely decays into the final 
states /i~A'g, fi~ K®, \i" p°, and pTu). Detailed GEANT simulations have been made for these decay modes. The 
calculated geometrical acceptances for the fi" p° and y~uj modes are very small. We therefore restrict ourselves to the 
H~Kg and y~Kl modes. We have not made an explicit event selection for the y~Kg and y~K\ modes. However, 
the GEANT simulations predict the fraction of these decays that would have passed our p — > /i~7 and p — > /i~7r° 
selection requirements. We can use these calculated efficiencies to set an upper limit on the fx~~ Kg and /i~ K^ modes. 
The calculated efficiencies e for these decays to satisfy the trigger and either the /i~7 or fi~Tr° search criteria described 
previously are (3.3 ± 1.0) x 10 -6 for the \i~ K s mode, and (5.0 ± 1.5) x 10~ 6 for the fj," K n L mode. Substituting these 
values into Eqs. 3 and 5 yields the limits Tp/B(p — > fi~K s ) > 4 x 10 3 yr and Tp/B(p — > y~ A"£) > 7 x 10 3 yr at 
90% C.L. To verify these results we have removed the last cut on our final kinematic distribution and used the more 
complicated statistical method described in Ref. |16[] . The resulting limits on the yT K s and y~ K Q L decay modes are 
consistent with our quoted results. 



V. SEARCH FOR ANTIPROTON DECAY WITH A FINAL STATE ELECTRON 

We begin with the sample of 37.8 x 10 6 events recorded with the Ei, a i ance , E to tai, and electron triggers described 



in Sec. [II B. Offline, after final calibration of the scintillation counters, a re-application of the upstream veto counter 
requirement reduces this sample to 34.2 x 10 6 events. The remaining backgrounds can be suppressed by requiring that 
the signals from the tracker, calorimeter, and scintillation counters are consistent with events arising from p — > e~ + X 
decays. In the following we first consider the simplest decay channel p — > e~7, and then describe the search for more 
complicated decay modes. 



A. p — > e 7 event selection 



To select events consistent with the decay p — ► e 7 we begin by requiring: 

(1) A pattern of hits in the dE/dx counters consistent with the presence of one, and only one, charged track (16.9 x 10 6 

events) . 

(2) A pulse height in excess of 1.5 MIPs in the preradiator counters (9.6 x 10 6 events). This is a first, very loose, 

electron candidate requirement. 

(3) No signal in excess of 0.6 MIPs in both a horizontal and vertical TC counter (7.6 x 10 6 events). Note that electrons 

and photons are expected to produce electromagnetic showers that are fully contained within the calorimeter 
and the TC lead wall, and do not penetrate into the TC. 

(4) Either 2 or 3 clusters reconstructed in the calorimeter (2.5 x 10 6 events). A fully contained p — > e~7 event would 

normally result in two calorimeter clusters. However a third cluster can arise if, for example, a large angle 
photon is radiated from the electron as it passes through the preradiator. 

(5) More than 90% of the total energy recorded in the calorimeter be outside of the cells immediately surrounding 

the beampipe (1.3 x 10 6 events). This requirement suppresses backgrounds associated with beam halo outside 
of, but close to, the beampipe. The requirement also ensures that a significant fraction of the energy associated 
with the observed electromagnetic showers has not "leaked" out of the calorimeter into the uninstrumented 
region within the beampipe. 

These very loose selection criteria were designed to reduce the sample to a reasonable size before running the track 
reconstruction code. After track reconstruction we impose the following requirements: 

(6) There must be one and only one reconstructed track that extrapolates to a calorimeter cluster (1.0 x 10 6 events). 

(7) The reconstructed track must be consistent with originating from within the fiducial volume of the decay tank. 

Specifically, the point of closest approach of the track to the beamline must be in the region < z < 375 cm 
(2.7 x 10 5 events). 

(8) The distance of closest approach of the reconstructed track to the beamline jl0| must be less than 1 cm (6.0 x 10 4 

events) . 
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(9) There must be a preradiator signal that exceeds 4 MIPs (2.2 x 10 4 events). This requirement was based on a 

study of the preradiator signals associated with the ir° — > 77 calibration data samples. 

(10) The electron candidate must be coplanar (Acf> > 3.05 radians) with the beamline and the highest energy neutral 
cluster (467 events). 

Figures [^a and^b compare, respectively, the transverse-momentum-balance and total-energy distributions for these 
467 events with predictions from a GEANT simulation of p — > decays. The observed distributions have peaks 
corresponding to large transverse-momentum imbalances and low total energies. In contrast to this, the simulated 
signal distributions peak at small transverse-momentum imbalance and at total energies corresponding to the beam 
energy. The backgrounds can therefore be further suppressed by requiring that candidate events are kinematically 
consistent with p — > e~7 decays. 



In general, a two-body decay in which the daughter particles have known masses can be completely specified by 
9 parameters (two momentum-vectors and a decay point xo,yo, zq). Consider p — ► e~j decays in which the parent 
antiprotons are traveling along the beamline with momenta po, energy Eq. In this case there are 6 constraints: (i) 
xo = 0, (ii) y Q = 0, (iii) E(e) + E(j) = E , (iv) p x (e)+p x (j) = 0, (v) p y (e) +p y (j) = 0, and (vi) p z (e) +p z {l) = Pa- 
tience, to completely specify ap-» e~j decay occurring within the APEX decay tank, we need to specify the values 
of 3 parameters, which we choose to be zq, E(e), and the ^-coordinate of the intersection of the electron trajectory 
with the calorimeter x ca i(e). In the following we will refer to a specific set of values for these parameters as a decay 
hypothesis. For a given decay hypothesis the intersections of the electron track with the 3 horizontal and 3 vertical 
tracking planes can be predicted, along with the positions and energies of the calorimeter clusters associated with the 
electron and photon. To quantify how well a given decay hypothesis describes the observed locations of the hits in 
the tracker planes, and the observed locations and energies of the clusters in the calorimeter, wc define: 



where the sum is over the 6 calorimeter measurements ( [x,y]-positions and energies for the electron and photon 
clusters) and (n — 6) tracking measurements. The A* are the differences between the measured and predicted values for 
the i th measurement, and the cr, are the corresponding measurement uncertainties. Note that the track reconstruction 
code requires that at least 4 of the 6 scintillating fiber planes record a hit along the track, and hence the number of 
measurements used in the x 2 sum is in the range 10 < n < 12. The best decay hypothesis for a given event is found 
by finding the set of parameters (zo, E(e), and x ca i{e)) that minimizes the \ 2 P er degree of freedom (% 2 /dof = X 2 /n). 
In this fitting procedure, hypotheses are rejected if the x 2 /dof > 10. Only 3 of the 467 events shown in Fig. ga have 
a best-fit decay hypothesis that survives this loose kinematic requirement. In Fig. ^a the distribution of the fitted 
values of x 2 /dof for these events are compared with the predicted GEANT Monte Carlo distribution for p — > e~j 
decays that pass the trigger and offline selection requirements. There is no evidence for a signal. The observed events 
have values of x 2 '/dof larger than would be expected for typical p — > decays. 

To set a limit on t/B we require that the x 2 /dof for the best-fit hypothesis is less than a maximum value, which 
we set equal to 2. No events survive this requirement. The efficiency for p — > e~7 decays to pass the trigger and offline 
selection requirements has been evaluated using the GEANT Monte Carlo program to generate 10 5 simulated decays 
uniformly along the beamline within the decay tank. We find e = (5.2 ± 1.4) x 10~ 4 , where the uncertainty on e arises 
from the systematic uncertainties on the trigger calibration, calorimeter energy scale, track reconstruction efficiency, 
and Monte Carlo statistics. The trigger and calorimeter scale uncertainties yield contributions to ay of ±14% and 
±21%, respectively, and have been evaluated by analyzing GEANT Monte Carlo p — > e~j samples in which the 
simulated trigger and calorimeter scales have been changed by ±1 a. The overall systematic uncertainty on e was 
calculated by adding these contributions in quadrature with the uncertainty on the track reconstruction efficiency 
(±7%) and Monte Carlo statistics (±1%). Inserting the calculated e and a e into Eqs. 3 and 5 we obtain the result 
Tp/B(p — > e~7) > 7 x 10 5 yr (90% C.L.). To verify this result we have removed the cut on the x 2 /dof distribution 
and used the more complicated statistical method described in Ref. |l(|. The resulting limit on Tp/B(p — > e~7) is 
consistent with our quoted result. 



B. p — > e 7 kinematics and a limit on t/B 




(6) 
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C. Search for p — > e 7r° 

To select events consist ent with the decay p — > e~ir° we begin by applying the p — > e~7 selection requirements 



(l)-(9) described in Sec. VA. These requirements select 2.2 x 10 4 events in which there is a candidate electron 
together with 1 or 2 additional clusters in the calorimeter. We would expect genuine p — > e~7r° decays to be either 
reconstructed as an electron plus one additional cluster (if the two photons from the tt° decay are not spatially resolved 
in the calorimeter), or as an electron plus two additional clusters (if the two photons are resolved). If there are two 
additional clusters, we require the cluster-cluster mass to be less than 500 MeV/c 2 (2.1 x 10 4 events). It is convenient 
to analyze all the p — ► e _ 7r° candidate events as two-body decays. To do this we reduce the electron plus two cluster 
events to an electron plus one cluster topology by combining the two neutral clusters. The combined cluster has an 
energy equal to the sum of the energies of the component clusters, and is located in the (x,y)-plane at a position given 
by the energy-weighted centroid of the component clusters. We now require that the electron candidate is coplanar 
(A<j> > 2.9 radians) with the beamline and the neutral cluster (1091 events). Applying the two-body kinematics 
fitting procedure described in Sec. VB, only 8 events have x 2 /dof < 10. In Fig. ^| the distribution of the fitted values 



of x 2 /dof for these events are compared with the predicted GEANT Monte Carlo distribution for p — > e~7r° decays 
that pass the trigger and offline selection requirements. There is no evidence for a signal. The observed events have 
values of x 2 /dof larger than would be expected for p — ► e~ir° decays. 

To set a limit on r / B we require that the x 2 /dof for the best-fit hypothesis is less than a maximum value, which we 
set equal to 2.5. No events survive this requirement. The efficiency for p — > e~7r° decays to pass the trigger and offline 
selection requirements has been evaluated using the GEANT Monte Carlo to generate 10 5 simulated decays uniformly 
along the beamline within the decay tank. We find e = (3.0 ± 0.7) x 10 -4 , where the uncertainty on e arises from 
the systematic uncertainties on the trigger calibration, calorimeter energy scale, track reconstruction efficiency, and 
Monte Carlo statistics. Inserting the calculated e and <r e into Eqs. 3 and 5 we obtain the result Tp/B(p — > e~7r°) > 
4 x 10 5 yr (90% C.L.). To verify this result, we have removed the cut on the x 2 /dof distribution, and used the more 
complicated statistical method described in Ref. Jig ]. The resulting limit on Tp/B(p — > e~ir°) is consistent with our 
quoted result. We have also searched for p — » e^ir decays by employing an analysis in which the two-body fitting 
procedure is replaced by a series of cuts on kinematic distributions [ff8| . This alternative analysis also yields a limit 
on Tp/B(p — » e~7r°) similar to our quoted result. 

D. Search for p — » e~r/, e~Kg, e~co, and e~jj 

We begin by describing the search for p — > e~r\. The selection requirements employed for this decay mode are also 
used fo r the s earch for the e~Kg, e~ui, and e~jj final states. The results for these latter decay modes are presented 



in Sec. VD2 



1. p — ► e rj 

The 77 decays into neutral modes (77 — > 27, 37r°, ir°2j) with a branching fraction of ~ 0.7. The decay p — > e~?7, with 
the rj subsequently decaying into these neutral modes, would give rise to an electron plus two-or-more photons in the 
final state. To select events consistent with the decay p — > e~r\ we begin with the p — > e~7 selection requirements 



(l)-(9) described in Sec. V A, which yield the sample of 2.2 x 10 events in which there is a candidate electron together 



with 1 or 2 additional clusters in the calorimeter. We then require that: 

(i) There are two clusters in addition to the electron candidate (1.2 x 10 4 events). 

(ii) The lowest energy neutral cluster has an energy exceeding 1 GeV (1222 events), and the other two clusters both 

have energies in the range 2 < E < 6 GeV (394 events). 

(hi) The electron candidate is coplanar (A</> > 2.8 radians) with the beamline and the neutral cluster pair (46 events). 
The coplanarity requirement is based upon a GEANT Monte Carlo study of simulated p — * e~ij decays within 
the APEX decay tank, and is less stringent than the corresponding requirements used for the p — > e~7 and 
p — > e _ 7T° searches described previously. 

(iv) The total energy measured in the calorimeter is in the range 7 < E tot < 13 GeV (30 events). 

(v) The transverse-momentum-balance variable p^ 1 is less than 0.3, where p^ 1 = \pt(e) — Pt{ll)\j '{pt(e) + Pt 

(24 events). 
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(vi) The event lies within the dashed-region of the {M x , Mi2)-plane shown in Fig. ^ (2 events), where the missing- 
mass-squared (M x ) is computed using the incoming p four-vector and the outgoing electron candidate four- 
vector, and Mi 2 is the measured mass of the cluster-pair. The region accepted contains 86% of the simulated 
p — ► e~rj decays that survive the previous selection requirements. 

The distribution of the observed events in the (M x , Mi2)-plane is shown in Fig. [7|a for the 24 events that survive 
the selection criteria up to (but not including) the (M x , Mi2)-requirement. The event population suggests that the 2 
events that survive the (M x , Mi2)-requirement are associated with the high-mass tail of the background distribution. 
However, we cannot distinguish these events from genuine p — e~r/ decays, and hence our limit on this decay mode 
is based on the observation of 2 events, yielding \x max — 5.7 at 90% C.L. The efficiency for p — > e~r\ decays to pass 
the trigger and offline selection requirements has been evaluated using the GEANT Monte Carlo to generate 6 x 10 4 
simulated decays uniformly along the beamline within the decay tank. We find e = (4.4 ± 1.1) x 10~ 5 , where the 
uncertainty on e arises from the systematic uncertainties on the trigger calibration, calorimeter energy scale, track 
reconstruction efficiency, and Monte Carlo statistics. Inserting the calculated e and a e into Eqs. 3 and 5 we obtain 
the result T ¥ /B(p -> e~r)) > 2 x 10 4 yr (90% C.L.). 



2. p — > e Kg, e u, and e 77 



The selection requirements defined for the p — > e~r\ search will also select candidates for other two-body p decay 
modes with an electron in the final state, provided the neutral particle accompanying the electron (i) is sufficiently 
massive to produce events that satisfy an (M x , Mi 2) -requirement, and (ii) decays at least some of the time into 
daughter particles that are observed as two neutral clusters. In particular, the p — > e~r] search results enable us to 
place limits on the decay modes p —* e~Kg, e~u>, and e~77. For these modes we require candidate events fall in the 
region defined by 0.1 < M x < 1.0 GeV 2 /c 4 and 0.3 < M 12 < 1.0 GeV/c 2 (Fig. 0). The corresponding limits will be 
based upon the observation of 4 events, yielding (i max = 8.9 at 90% C.L. Note that, of the simulated signal events 
that pass all the prior selection requirements, the final (M x , Mi 2 )-requirement accepts 45% of the simulated e~Kg 
decays, 70% of the e~u decays, and 57% of the p — > e~77 decays. The efficiencies for these decays to pass the trigger 
and offline selection requirements have been evaluated using the GEANT Monte Carlo program, and are given by 
e = (3.0 ± 0.9) x 10" 6 for the e~ K° s mode, e = (6.8 ± 1.9) x 10" 7 for the e~u mode, and e = (4.9 ± 1.3) x 10~ 5 for 
the e~77 mode. Inserting the calculated values for e and a e into Eqs. 3 and 5, at 90% C.L. we obtain the results 
Tp/B(p — * er X) > 900 yr for the e~Kg mode, 200 yr for the e~u> mode, and 2 x 10 4 yr for the e~77 mode. 



E. Search for 



A K® L with an energy of a few GeV will travel of order 100 m before decaying. Hence, p — » e~ K\ decays occurring 
within the APEX fiducial volume would be expected to produce an electron accompanied by one neutral cluster in the 
calorimeter. In our search for p — > e~Kl decays we begin with the p — ► e~7 selection requirements (l)-(9) described 



in Sec. VA, which yield the sample of 2.2 x 10 events in which there is a candidate electron together with 1 or 2 



additional clusters in the calorimeter. We then require that: 

(i) There is only one cluster in addition to the electron candidate (9300 events). 

(ii) The electron candidate is coplanar (A<f> > 3 radians) with the beamline and the neutral cluster (300 events). 

(iii) The total energy measured in the calorimeter is in the range 7 < E to t < 13 GeV (157 events). 

(iv) The electron cluster and the neutral cluster both have energies in the range 2 < E < 6 GeV (11 events). 

(v) The missing-mass-squared lies in the range 0.1 < M x < 0.4 GeV 2 /c 4 , where M x is computed using the incoming 

p four-vector and the outgoing electron candidate four-vector (3 events). 

The Mj^-distribution for the observed events is shown in Fig. || for the 11 events that survive the selection criteria 
up to (but not including) the M x requirement. The event population suggests that the 3 events that survive the M x 
requirement are associated with the high mass tail of the background distribution. However, we cannot distinguish 
these events from genuine p — > e~ decays, and hence our limit on this decay mode is based on the observation 
of 3 events, yielding [i max = 7.3 at 90% C.L. The efficiency for p — > e~K*[ decays to pass the trigger and offline 
selection requirements has been evaluated using the GEANT Monte Carlo to generate 10 5 simulated decays uniformly 
along the beamline within the decay tank. We find e = (2.2 ± 0.5) x 10 -5 , where the uncertainty on e arises from 
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the systematic uncertainties on the trigger calibration, calorimeter energy scale, track reconstruction efficiency, and 
Monte Carlo statistics. Inserting the calculated e and a t into Eqs. 3 and 5 we obtain the result Tp/B(p —> e~ > 
9 x 10 3 yr (90% C.L.). 

VI. CONCLUSIONS 

A search has been made for 11 two-body and 2 three-body decay modes of the antiproton in which there is an 
electron or muon in the final state. No statistically significant signal indicating p decay has been observed. Our limits 
on t/B are summarized in Tab. |[ and are significantly more stringent than the limits currently quoted by the Particle 
Data Group In particular, our results place the first explicit experimental limits on the muonic decay modes of 
the antiproton, and the first limits on the decay modes e _ 77 and e~uj. Our most stringent limit is on the decay mode 
p — > e _ 7 for which we find r/ B > 7 x 10 5 yr (90% C.L.), which is an improvement by a factor of 400 over the prior 
T861 result. Noting that our limits for the simplest two-body decay modes are in the range 10 4 to 10 5 years, Eq. |l| 
suggests that if there is a CPT violating interaction that gives rise to these decay modes, and that is described by a 
dimension-6 operator, the associated mass scale uij is greater than O(10 9 ) GeV/c 2 . This is well beyond the reach of 
direct searches at high energy colliders. 
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TABLE I. Summary of results: 90% C.L. limits on t/B for 13 antiproton decay modes. 



Decay Mode 


t 1 d ijimit 


Decay Mode 


t j d ijimit 




(years) 




(years) 


/x~7 


5 x 10 4 


e~7 


7 x 10 b 




5 x 10 4 


e -7T° 


4 x 10 5 


fi-r] 


8 x 10 3 




2 x 10 4 


^-77 


2 x 10 4 


e _ 77 


2 x 10 4 




7 x 10 3 


e-Kl 


9 x 10 3 




4 x 10 3 


e-K% 


9 x 10 2 






e u) 


2 x 10 2 
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FIG. 1. Schematic of the APEX detector. 
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FIG. 2. The four calorimeter quadrants used for triggering. 
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FIG. 3. Distribution of neutral cluster energies for the 14 events that pass the p — > /n 7 selection criteria described in the text 
(open histogram) compared with the predicted distribution for a signal corresponding to t^/B — 5000 yr (hatched histogram). 
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FIG. 4. Distribution of event masses for the 16 events that pass the selection criteria described in the text (open histograms) 
compared with predictions from a GEANT simulation (hatched histograms) for the decays (a) p — > /u - 7r°, (b) p — > fi^r/, and 
(c) p — ¥ /i _ 77- The two entries in the lowest bin are events with calorimeter energy less than the -k rest energy. The predicted 
signal distributions are normalized to correspond to t^/B = 5000 yr for the fi~iv mode, 1000 yr for the [i~r] mode, and 2000 yr 
for the /i - 77 mode. 
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FIG. 5. Kinematic distributions for the 467 p —> e~y candidate events prior to the kinematic analysis described in the text. 
The observed distributions are compared with GEANT Monte Carlo predictions for p — > e _ 7 decays normalized to correspond 
to a lifetime of 5000 years, and shown for (a) a measure of the transverse momentum balance in the event, where -Pt(I) and 
Pt(2) are the transverse momenta of the two highest energy clusters in the event, and (b) the sum of the energies of the two 
highest energy clusters in the event. 
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FIG. 6. Distribution of \ 2 P er degree of freedom for (a) the 3 p -> e~7 candidate events surviving the kinematic fitting 
described in the text, compared with the GEANT Monte Carlo predictions for p —> e~7 decays normalized to correspond 
to a lifetime of 100,000 years, and (b) the 8 p — > e~7r° candidate events surviving the kinematic fitting described in the 
text, compared with the GEANT Monte Carlo predictions for p — ► e~7r° decays normalized to correspond to a lifetime of 
100,000 years. 
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FIG. 7. Observed event populations in the (a) missing-mass-squared versus cluster-cluster mass plane for candidate events, 
and (b) simulated p — > e~r\ decays that pass the prior selection requirements described in the text. The simulated data 
correspond to a lifetime of 600 years. The dashed lines delineate the region accepted by the analysis cuts for p — > e~7/ decays 
and the solid lines delineate the region accepted by the analysis cuts for other modes, as described in the text. 
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FIG. 8. Distribution of missing-mass-squared for the 11 p — > e~ K*l candidate events that survive all of the requirements 
described in the text except for the missing-mass-squared requirement. The observed distributions are compared with the 
GEANT Monte Carlo predictions for p — ► e~ ' K% decays normalized to correspond to a lifetime of 2,500 years. 
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